Strong electron-photon coupling in one-dimensional quantum dot chain: 

Rabi waves and Rabi wavepackets 
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We predict and theoretically investigate the new coherent effect of nonlinear quantum optics - 
spatial propagation of Rabi oscillations (Rabi waves) in one-dimensional quantum dot (QD) chain. 
QD-chain is modeled by the set of two-level quantum systems with tunnel coupling between neigh- 
boring QDs. The space propagation of Rabi waves in the form of traveling waves and wave packets is 
considered. It is shown, that traveling Rabi waves are quantum states of QD-chain dressed by radi- 
ation. The dispersion characteristics of traveling Rabi waves are investigated and their dependence 
on average number of photons in wave is demonstrated. The propagation of Rabi wave packets 
is accompanied by the transfer of the inversion and quantum correlations along the QD-chain and 
by the transformation of quantum light statistics. The conditions of experimental observability are 
analyzed. The effect can find practical use in quantum computing and quantum informatics. 
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I. INTRODUCTION 



Rabi oscillations are periodical transitions of a two- 
level quantum system between its stationary states un- 
der the action of an oscillatory driving field, see e.gr^. 
The phenomenon was theoretically predicted by Rabi 
on nuclear spins in radio-frequency magnetic field 3 - and 
was firstly observed by TorrejA Afterwards, Rabi os- 
cillations were discovered in various physical systems, 
such as electromagnetically driven atoms 5 - (including the 
case of Rydberg atomic states 6 -), semiconductor quan- 
tum dots (QDs) 7 , and different types of solid-state 
qubits (superconducting charge qubits based on Joseph- 
son junctions&2Jii, spin-qubitsii, semiconductor charge- 
qubits 12 ). In real physical systems the ideal picture 1 
of Rabi effect can be essentially modified by additional 
features, such as the time-domain modulation of the 
field-matter coupling constan t 13 ' 14 , the phonon-induced 
dephasing 15 and the local-field effect o 17 ' 18 ' 19 . New phe- 
nomena appear in Rabi oscillators with broken inver- 
sion symmetry 2 ^ and in systems of two coupled Rabi 
oscillator a 21 i 22 i 23 i 24 . 25 i 26 i 27 . 

In spatially extensive samples with a large number of 
oscillators the propagation effects come into play. As a 
result, the mechanism responsible for Rabi oscillations 
causes also a number of nonstationary coherent optical 
phenomena, such as optical nutation, photon echo, self- 
induced transparency, eto 2 ^. In low-dimensional systems 
the propagation effects also take place. For example, the 
numerical modeling of the coherent intersubband Rabi 
oscillations in a sample comprising 80 AlGaAs/GaAs 
quantum wells^ shows that the population dynamics de- 
pends on the quantum well position in the series. This 
result demonstrates strong radiative coupling between 
wells and, more generally, significant difference in the 
picture of Rabi effect for single and multiple oscillators. 
Another aspects of such a difference, namely, the effects 



of quantum interference and correlations between pho- 
tons in multiatom fluorescence, are demonstrated i n 30 i 31 
(on the example of atomic chains in nanofibers). From 
practical point of view, the effect of Rabi oscillations is a 
key ingredient for realization of binary logic and optical 
control in quantum informatics and quantum computing. 

The theoretical analysis of Rabi oscillations is highly 
diversified both in form and content. The common fea- 
ture is the impossibility of consideration of electromag- 
netic field influence as a small perturbation. Different 
ways of description are used in the analysis of the model 
problem of Rabi oscillations in a single two-level atom 1 . 
One such a way, the probability amplitude method, con- 
sists in solving of the Schrodinger equation for wave- 
function |\&), which is the superposition of various atom- 
photon states. The second one is the Heisenberg operator 
method which is based on the analysis of the photonic 
and atomic operators time evolution. And finally, the 
third way is the unitary time-evolution operator tech- 
nique. As is demonstrated ini, all these methods lead 
to identical solutions; the choice of the concrete one is 
determined by the convenience considerations. 

Taking into account the quantum nature of the electro- 
magnetic field has the special significance in the analysis 
of Rabi oscillations in complex systems. Two different 
cases can be marked out. The first could be called qua- 
siclassical one 2 ^. In this case external field has classi- 
cal nature. However, the true field is supposed to be 
concordant with the quantum motions of the particles. 
So, the field should contain the contribution of the in- 
duced polarization, which has a quantum nature. An- 
other case is really quantum and takes into consideration 
the photon structure of the electromagnetic field 1 . Dur- 
ing the process of energy-level transitions in the atom the 
photon structure of the field is also transformed, there- 
fore the atom-photon dynamics should be considered self- 
consistently. 

In this paper we build for the first time a theoretical 
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model of a distributed system of coupled Rabi oscillators 
and predict a new physical effect: propagation of Rabi 
oscillations in space in the form of traveling waves and 
wave packets. Quantum oscillators in the system that 
exhibit Rabi-wave propagation interact strongly enough 
to stare and exchange the e-h pair excitation between 
light emission and light absorption acts. Often in quan- 
tum optics an "all-matter" picture is employed, where 
the dynamics of electromagnetic field is integrated out, 
for example in the optical Bloch equation a 32 ' 33 . I n 34 i 35 
the light-matter interaction is treated in an "all-light" 
picture (Lippmann-Schwinger equation approach). 

In our model the quantum nature of light is also fully 
accounted. As the theoretical approach the probability 
amplitude method, generalized for the case of ID-chain, 
is convenient. 

The paper is organized as follows. In Sec|TT]the model 
is formulated and equations of electron-photon dynam- 
ics both for discrete chain and continuous limit are ob- 
tained. In Sec lIIII the plane Rabi waves are considered 
(their dispersion equations, eigenmode structure, disper- 
sion curves properties). Sec[lV] is devoted to the Rabi 
wavepackets investigation, their spatial-temporal dynam- 
ics is analyzed, experimental observability of Rabi waves 
is discussed. In Sec|V]the Rabi waves in the case of classic 
light are examined and distinctive features of scmiclassi- 
cal consideration are treated of. In Scc lVII thc role of the 
local-field effects and conditions of this effects negligibil- 
ity are discussed. In Sec lVIII the space-time structure of 
electron-electron and electron-photon correlators is con- 
sidered. The main results of the work are formulated in 
SecEnH 

II. MODEL 
A. Hamiltonian 




FIG. 1: Schematic picture of the QD-chain interacting with 
single-mode electromagnetic field (a) and the energy levels 
diagram for p-th QD (b). 

Consider an infinite periodical one-dimensional chain 
of identical QDs, containing an electron in the size- 
quantized conduction band. It is assumed that in each 
p-th QD there are at least two one-electron orbital 
states, ground \b p ) and excited \a p ), with transition fre- 
quency loq. Neighboring QDs are coupled via electron 



tunneling 27 , i.e. the electron can go from state \a p ) 
to state |a p ±i) and from \b p ) to state \b p ±\). Transi- 
tions between ground state and excited state belonging 
to different QDs are neglected: (a p ±i \b p ) w 0. Let 
the QD chain be exposed to a plane wave of quantum 
light, which is incident on the chain at an angle a, 
E(x) = £(ae lkx + a + e~ lkx ), where a, a + are photon cre- 
ation and annihilation operators, £ = /Vbe, Vq is 
the normalizing volume, e is the unit polarization vector, 
wavenumber k = wcosa/c. The dependence of the field 
on transverse coordinates is negligibly small because of 
electrical smallness of QDs. The case of oblique propa- 
gation is of interest due to the possibility of wavenumber 
variation with variation of incidence angle a at constant 
frequency to. Then the Hamiltonian of the system in the 
rotating wave approximation 1 reads 

H = H d + H f + H df +H T + AH, (1) 

where H d = (Hujq/2) J^p^zp i s Hamiltonian of QD-chain 
in the absence of electron tunneling and QD-field inter- 
action, a zp — \a p ) (a p \ — \b p ) (b p \, Hf = huia + a is the 
Hamiltonian of the free electromagnetic field 3 ^. The com- 
ponent 

H df = hg J2 (t+ae lk P a + a p a+ e~^ a ) (2) 
p 

describes QD-field interaction, where g — —iiE/h is in- 
teraction constant, fi is the QD dipole moment. Note 
that all transition dipole moments in the chain are as- 
sumed to be real values and to have fixed orientation. 
The operators a p = \a p ) (b p \ and a~ = \b p ) (a p \ are the 
transition operators from ground to excited states and 
vice versa for p-th QD. The term 

H T = -ttiii i\a P ) (ap+il + K) (%>-i|) 

(3) 

p 

takes into account interdot electron tunneling; £1.2 are 
the electron tunneling frequencies for the excited (£1 ) and 
ground (£2) states of the QDs. The term AH corresponds 
to the local-field effects originated from the dipole-dipole 
electron-hole intradot interactio n 17 ! 18 ' 19 ; in the mean- 
field approximation this term is given by 

47T 

AH = V »(Np,) ]T(<V(a, + ) + *+<*-)). (4) 
p 

Here N_ is the depolarization tensor (see Eq. (18) in 
Refjl 9 .), V is the volume of QD. 

B. Equations of motion 

The state vector of the " QD-chain+light" system may 
be represented in terms of the eigenstates of isolated QDs 
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Here \b p ,n) — \b p ) ® |n), \a p ,n) — \a p ) (8 |n), where 
|n) is the light Fock state with n photons, i? Pi „, A p , n 
are the probability amplitudes. In the interaction pic- 
ture the evolution of the system is described by the 
nonstationary Schrodinger equation ihd t \^) — V\^) : 
where the interaction Hamiltonian has the form V — 
exp(iH f t/h)(H d + H df + H T + AH)exp(-iH f t/h). This 
Schrodinger equation leads to the following equations for 
the probability amplitudes 



dA Pi n iujQ . 

— 7^ — — ^-A p ^ n + lt;i (A p _i^ n + A p+ \, : 



(0) 



m 

dBp, n +l _ iuJ R A. it f R I d ^ 
— — £>p,n+l + \£> P -l,n+l + £> P +l,n+l) 

- igV^+lA p . n e- l ( kpa -^ - iAuiA Pin+1 J2K,m B P,m, 



where 



4tt ~ 
Alu = j^ViKv) 



(8) 



is the local- field induced depolarization shift 19 . Obtain- 
ing this equations we have taken into account that the in- 
teraction @ can cause the transitions between the states 
\a p ,n), |b p ,n+l) only. As it is seen from Eqs.((6])-([7]), 
two competitive mechanisms manifest themselves in the 
light - QD-chain coupling: the local-field induced nonlin- 
earity and the dispersion spreading due to the tunneling. 

To this point we have taken no account of the decay 
processes inside the QDs. To take this processes into ac- 
count one should allow for processes of interaction of e-h 
pair in QD with phonon bath. Such interaction can be 
described by conception of quantum trajectorie o 37 ' 38 ! 39 
of e-h pair, which are the superpositions of deterministic 
evolutions and random jumps under the action of Lind- 
blad operators. In this case the Hamiltonian ([I]) should 
be replaced by the efficient non-Hermitian Hamiltonian 
H e ff = H — ihJ2 p X p X p , where X p is Lindblad operator 
for p-th QD. Suppose that the density of phonon states of 
the bath has the form of Lorentz line with frequency u>q 
and width A. If the bath is weakly coupled with QD then 

X p = ^\/2a p (see Refs 3 -^^). It means that the decay 
can be taken into account by substitution luq ujq — iX 
into Eq[B] and ujq — > ujq + iX into Eq[71 Let us restrict 
here consideration to linear regime of the carrier motion 
and omit the terms O(Aoj) in ©-([I]). In the Sec I VII we 
shall analyze the limitations imposed by omission of non- 
linear terms. By this means, omitting the terms 0(Auj) 1 
we arrive at the set of difference-differential equations for 



vectors \P p , n (t) 



B 



Ap : n(t) 
(*) 



as follows: 



2 o z - -1 - ig\pn^\k p (t) 



(9) 



-t|(*p_l, n + *p+l,„) : 



where / denotes 2D unit operator, k p (t) = a x exp[i(u;i — 
kpa)a z ], | = [(6 + + (£i - 6)<7»]/2. Coefficients in 
([9|) are expressed in terms of Pauli matrices 
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acting on vectors *& P:n (t). 



C. Continuous limit 

Equations © are the recurrent ordinary differential 
equations. In some cases it is more convenient to trans- 
form them into the system of partial differential equa- 
tions. To do this we should turn to the continuous limit, 
making the standard substitutions (see, for example^) 

a 2 d 2 ^f n , h p — > k(x). Then the system ([9]) leads to 



2^ 



ILOq 



X 



I — igyn + 1 k(u)t — kx) 



-ia 2 id^ r 



(11) 

Formula (jlip represents the system of partial differen- 
tial equations with the variable coefficient k — k(ujt~kx). 
This is the basic system for most of our further calcula- 
tions. 



III. TRAVELING RABI WAVES 

A. Eigenmodes 

Consider elementary solution of system 
(jll[) in the form of damping traveling waves: 

A n = Une i(h+k/2)x e -i(v+u/2)t e -\t/2 t Bn+i = 

v n j rl e %(h ~ k l' r)x e~ il - v ~ u 'l T)t e~ xt l' 1 , where ft, is a given 
wave number and v is the eigenfrequency to be found, 
u n , v n are the unknown constant coefficients. Substi- 
tuting them into equations (jlip and omitting common 
factors, we obtain 



{v - A/2 + i?i (h)} u n ~ gVn + Tv n+1 =0, (12) 



gVn + lu n - {v + A/2 + & 2 {h)} v n+1 = 0, (13) 

where A = luq — to, 

$ia{h)=^ 2 [2-a 2 {h±k/2)\ (14) 

System (TT2")) ~ (ri"3l has a nontrivial solution if its deter- 
minant is equal to zero, i.e. we obtain quadratic equation 
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which connects v and h. Solving the equation with re- 
spect to v, we determine the eigenfrequencies of system 
as 

is 12 (n,h) = -^[Mh) + Mh)Tn n {h)}. (15) 

Here 



n n {h) = ^a^ + v^ + i), (i6) 

A eff {h) = A-# 1 (h) + Mh)- (17) 



It follows from (fT5|) that there are two eigenmodes for 
each value of photon number n, namely 

\^i,n(t)) = ^2 (At n (pa, t) \a p , n) + B i<n (pa, t) \b p , n + l)), 
p 

where i = 1, 2 is the number of the eigenmode, and coeffi- 
cients Ai -n (x, t), Bi tn (x, t) are defined by the expressions 



A 1<n (x,t) 



V X - A/2 + 7?! 

B hn+X {x, t) = c l e i{ - h - k W x e- i{ - Vl - u l 2 - i W t (19) 
and 

A 2 , n (x, t) = c 2 e l(h+k/2)x e- l{ " 2+u)/2 - lX/2)t , (20) 



B 2 ,n+l(x,t) 



C 2 g^/n + \ e i (h~k/2)x e -i{v 2 -u /2-i\/2)t 

v 2 + A/2 + i? 2 ' 



where C\ j2 = Ci t2 (n, h) are normalizing constants, i?i 2 = 
$1,2(^)1 the values V\ 2 = v\ 2 {n,h) are expressed by the 
Eqs. ([HI), HU), respectively. Each mode corresponds to 
the electron-photon entangled state, which partial ampli- 
tudes oscillate both in time and space. The expression 
(fTS]) with (fT9| or (|20| describes the siaie 0/ radiation- 
dressed QD-chain since the entangled state includes field 
states with different photon numbers (n and n+l). This 
states are the generalization of the single atom dressed 
state: 2 for the case of distributed system. The qualita- 
tive distinction of this case is the space-time modulation 
of dressing parameter: it propagates along the QD-chain 
according to the traveling wave law exp[j(fc:r — u)t)]. En- 
tanglement and dressing are caused by the interaction of 
light with QD chain and vanish in the limit of g — > 0. In 
that case, (fTTj) describes QD-chain electron-hole pairs in 
equilibrium and inverse states, respectively. 

Spatial oscillations of the partial amplitudes are due to 
QD-coupling. They vanish in the limit of £1.2 — » 0. Gen- 
erally, both the modes are excited simultaneously. How- 
ever, any of them can be excited separately by a proper 
choice of initial conditions. Similarly to the other coher- 
ent excitations in condensed matter, Rabi waves (|19p and 
(|20p introduce a new family of quasi-particles (we name 
them rabitons). One can apply to them the standard sec- 
ondary quantization technique. Similarly to the effect of 



the self-induced transparency, the Rabi wave propagation 
can be interpreted as the motion of a precessing pseudo 
dipole^. However, the coherence mechanisms in these 
two cases are principally different: in the Rabi wave the 
coherence is settled by the dispersion relation (fT5|) while 
in the case of the self-induced transparency it has soli- 
tonic character. 

Note that the eigenmodes (fT9|) and ([20]) each comprise 
traveling waves with different wave numbers h ± k/2. 
Physically, it means that the Rabi wave propagates in 
an effective periodically inhomogeneous medium formed 
by spatially oscillating (with period 2ir/k) electric field. 
Therefore, the diffraction appears in the system. In the 
limit k — ► 0, the medium turns homogeneous and the 
diffraction effect vanishes. Reflections of Rabi waves and 
their mutual transformations at the field inhomogencitics 
become possible. Thus one obtain a unique possibility to 
control the processes of the reflection and transmission of 
Rabi waves by varying the spatial structure of the light . 



B. Dispersion characteristics 

Note, that for the waves under consideration the reci- 
procity conditions break down: Vx, 2 {n,h) 7^ Vx t2 {n,—h). 
This is due to the presence of preferential direction, which 
is determined by the direction of photon mode propaga- 
tion along the QD-chain (sign of the value k). The form 
of the dispersion characteristics for typical values of pa- 
rameters is presented on Figf2] To interpret them one 
should note that in the limit g — > both of branches 
have cross points, two in general case and one at £1 = £2- 
Modes under examination are characterized by continu- 
ous spectrum (the value h varies continuously). 

The values u\ >2 are real for all real h. Therefore, the 
system being investigated is stable^-. Dispersion equa- 
tion for (rr^ - lTTn)) can be solved with respect to wave num- 
ber h(y) for given real-valued frequency v. In this case 
for every mode there are such critical values vf 2 that for 
v < vl r 2 the values h\^ 2 (v) become complex. Physically, 
it means opacity of the QD-chain for Rabi waves with 



such frequencies 41 . For v 2 



< v < v\ 



the QD-chain is 



opaque only for one of the modes and for v < v1 r 2 it is 
non-transparent for both of them. One can choose the 
system parameters in such a way that in some range of 
v the all real solutions h{y) are either negative or pos- 
itive (as an example see Fig[2ji). This implies one-way 
opacity of the system (Rabi waves could propagate in one 
direction only), which is caused by nonreciprocity of the 
QD-chain. The critical frequencies are determined from 
the condition vf r 2 (n) — v\, 2 {n, h^), where 



dvi. 2 (n, h) 



dh 



(21) 



h=h\ 



It is significant, that critical values vf r 2 depend on pho- 
ton number n (see Fig|3]). This gives an important result, 
namely: the conditions of transparency of the QD-chain 
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FIG. 2: Dispersion curves for Rabi waves for different param- 
eters values, a) £1 = £2 = 8g, A = 3g, ka = 1, n = 5; b) 
£1 = Wg, £2 = 3g, A = 0, ka = 0.143, n = 2; c) £1 = IO5, 
£2 = 0, A = — 5<?, fea = 1, n = 5. 



are different for Rabi waves with different photon num- 
bers. Note, that vf increases and z/T decreases with 
increasing of n. 

Usually both Rabi waves demonstrate the normal dis- 
persion and propagate in the same direction. But for a 
special choice of the parameters the abnormal dispersion 
region for one of the modes appears. For this mode phase 
and group velocities are oppositely directed (see FigJ^.). 
The dispersion characteristics have special behavior for 
absolutely nontransparent ground-state barrier: £2 = 
(Figf2fc). In this case the group velocity v gr = dv/dh of 
one of the waves is very small for all h, except for the 
narrow vicinities of the cross-points. Due to this fact 
the modes exchange their places when passing the cross 
points (e.g. on Figj^t v gr w for i^-mode between cross 
points —2 < ha < 1 and for ^2-mode at all other values 
of h). Thus, in current case for a wide range of h one of 



FIG. 3: Dependence of critical frequencies of Rabi waves on 
the photon number, for the input parameters as follows: £1 = 
10ff, £2 = 7g, A = 2(£i - |a) + &a 2 k 2 , ka = 1. 



the modes practically does not transfer the energy along 
the QD-chain. 



IV. RABI WAVEPACKETS 
A. Electron-photon dynamics 

To find the general solution of the system ([IT]) 

let us introduce the new variable $„(x,t) = 
e i{w t-kx)d z /2 e \t/2^ n ^ x ^ t y It leadg tQ the system of par- 
tial differential equations with coefficients independent of 
QD-position x: 



d t $n - i 



(2 - a 2 fc 2 /4) £ - gVn~+lx(t) 



+a A ka z td x & n - ia^d^ n = 0, 



(22) 



where x(t) = a x exp(— ia z At). This system can be solved 
by using the Fourier transform with respect to x: 



(23) 



Substituting Eg. (f2"3")l into Eq. (|2"2"|) . we obtain the system 
of ordinary differential equations 



-* n = t(£- 5 V£+Tx(t))* n , 



(24) 



where 3 = d(h) = [[#i(h) + ®i{h))I + (tfi(/i) - 
'&2{h))a z ]/2. The system (pM)) can be easily integrated, 
the solution is given by 



s n (M) =p„(M)* n (M) > 



where 



p n (h,t) = e il^i(h)+^(h)]t e i» z At e i[m&n n (h)]t^ 



(25) 



(26) 
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& = ipxi^y,^z)i m is the unit vector of the form 
(—2gy / n + l/O n , 0, A e ///fi„). The expression (|26)l can 
be written in matrix form as follows: 



p n (h,t) 



(f.-(h,t)e iS +^ t ip n {h,t)e i5 +^ t 
ip n (h, t)e lS - 93+ (h, t)e tS - W 



<Pn(h,t) = cos 



Q n (h)t A eff (h) ^ Q n (h)t 



Q n (h) 



sin ■ 



(27) 



(28) 



, 2gVn + l . n n (h)t 

Mh > t) = - i ^uKr sm ^^> (29) 

<5±W = ^[^iW + ^W±A], (30) 

Vt n (h) and A e ff(h) are determined by (fH)|) and (jT7)l . cor- 
respondingly, and <fr(/i, 0) can be obtained from initial 
conditions by the inverse Fourier transform: 



*«(M) = 7^ / e-« hI+a > k ^ x * n (x,0)dx. (31) 

Z7T 



Thus, for the vector of probability amplitudes *& n (x, t) 
we have the following expression: 

OO 



Expression (j32|) is the central result of this section. It 
allows one to characterize such observables as inversion, 
as well as the exciton-exciton and exciton-photon corre- 
lators for different initial states of light and QDs. 



B. Dispersionless approximation 

In this section we will analyze the propagation of Rabi 
waves with spatially localized initial state. The QD- 
chain issupposed to be initially either in the ground 
state (A n (x,Q) = for all n), or in the excited state 



(B n (x,0) = for all n), or in the mixed state (arbitrary 
superposition of ground and excited states). 

We describe the spatial-temporal dynamics of Rabi os- 
cillations by the spatial density of the inversion (the in- 
version per single QD) 

w(x,t) = aJ][K(x,i)| 2 - \B n+1 (x,t)\ 2 }. (33) 



Space-averaged temporal evolution of Rabi oscillations is 
characterized by the integral inversion 



w{t) 



w(x, t)dx. 



(34) 



Let us approximate the initial spatial distributions 
by the Gaussian beams Ca,b exp[— (x — cIa,b) 2 /2c^ b ], 
where Ca,b are normalization constants, cIa,b, &a.b are 
the positions of the beams and their widths, respectively 
(the indexes A, B refer to excited and ground state of the 
electron in the QD-chain, respectively). 

The Rabi wavepacket evolution can be qualitatively 
illustrated with the help of the dispersionless approxi- 
mation for Eqs. (|3"2")l . According to (|31[) we obtain for 
a„(h) 



#„(fc,0) = 



'M+l 



(h) 



a n (h) ~ ^ n)e -(h+k/2fal/2 e - ldA (h+k/2)^ (35) 



b n+1 (h) ~ C(n+ l) e -(h-k/2fa%/2 e ^d B (h-k/2)^ (gg) 

where c(n) is an arbitrary photonic distribution. One can 
see from (l3"5T) - (|3l>|) that the spatial spectrum of ^ n (x,0) 
is given by two Gaussian peaks with peak widths 1/cta.b 
and positions of the peak centers — ±fc/2, respec- 
tively. If a a,b are sufficiently large, the largest contri- 
bution to the integration in l|32p comes from the narrow 

vicinities of points h^\. Assuming di(h) = , di{h < f\) + 

(dtii/dh)\ h=h m (h- h[°l) and p n (h,t) = p n (h^l,t) (the 

dispersion effects neglecting) and doing standard calcu- 
lations, we obtain from Eg. ([32)1 : 



A n (x, t) = e ilkx-u>t)/2 e -\t I ]^ An ^ x + v + t ^ ^+ e ie+t + A ^ x + v - tj o)f- e ifl > * 
+ [B^ix + v+t^O^e' 9 ^ +B n+1 (x + v^t,0)f ] +e t ^ t ] e* kx ' 2 } , 
B n+1 (x, t) = e -H**-vt)/2 e -\t I An ( x + v + t ^ o^-e^i* + A n (x + v-t, O^+e^"* 
+ [i? n+1 ( a ; + ^,0)Cre<* + i?„ + i(x + U rt,0)C 1 + e^*] e lkx ' 2 } , 



ikx/2 



-ikx/2 



(37) 



(38) 



where the velocities v\ 2 are defined as 



vf = v ± (h$ >) ) = -Zid 4 k 



n n (fc/2)TA e// (fc/2) 

n„(fc/2) 



v ± „± f/l <oK 2 M-k/2)±A eff (-k/2) 
v 2 -v (h 2 )-&o * 0n (-fc/2) 



(40) 



(39) 
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and 9f 2 and Ci^i r ii2 are introduced to denote the fre- 
quency shifts and amplitude factors, respectively: 



1 



,(») 



(0)m _.,±/ fc (0)> 



_ n n (ftW)±A ef/ (ft^) 

^ 20„(/4° 2 ) 



(41) 
(42) 

(43) 



As is apparent from ([3"T)) - ([3"8)) . generally, any probabil- 
ity amplitude in the Rabi-wave packet is the superpo- 
sition of four subpackets. Two of them (first and sec- 
ond terms in expressions ([57)) - ([55)) ) correspond to the 
excited initial state and two another (third and fourth 
terms in ([37)1 - ([38)) ) correspond to the ground initial state. 
The partial subpackets in Eqs. ([37)) - ([38)) are character- 
ized by the different frequency shifts 9f 2 an d different 
velocities of motion vf 2 , given by ([39)) , ([40]) . This ve- 
locities coincide with the group velocities of traveling 
Rabi waves (TT9")) - ([2"0")) when h = ±fe/2: one can ele- 
mentary verify, that vf = dvi.2/ dh\ h=k ^ 2 = wj ,r 2 (fc/2), 



'h = dv ifl/ dh \ h =- 



k/2 



gr 



(—k/2). It is essential that 



velocities vf 2 as wei l as the frequency shifts 9f 2 depend 
on photon number n. It means that the spatial prop- 
agation of wavepacket is accompanied by the change of 
quantum light statistics (for example, in initially coher- 
ent light incoherent component appears). 

Under the some specific conditions the number of sub- 
packets could decrease. Two mechanisms of such de- 
crease are possible: the first one is the tending to zero the 
subpacket amplitude and the second one is the confluence 
of the subpackets velocities. 
As is seen from (t351>— (f3Dl). if 



(44) 



A eff (h?h = 0, 



for one pair of the subpackets the velocity synchronism 



condition 



is fulfilled. Notice that Eq. (gH is 



the analogue of the ordinary synchronism condition in 
single two- level system A = 0. The velocity synchronism 
condition can be fulfilled by fitting the value of the de- 
tuning A. If this is the case, there are three subpackets 
instead of four them due to the confluence mechanism 
acting. Note, that A e //(fc/2) and A e //(— k/2) can not 
be equal to zero together at k 7^ 0. 

If the system is initially prepared in the stationary 
state, and the synchronism condition is fulfilled in the 
corresponding point (for example B n+ i(x,0) = for any 
n and A e //(— k/2) = 0), only one subpacket is preserved. 
In this case the expression for the inversion density be- 
comes rather simple: 

w(x,t) = a^2A 2 Jx + £ 2 a 2 fci,0)[l - 2 sin 2 (gVn + It)}. 

n 

(45) 



Expression (|45|) sums up contributions of different pho- 
tonic states (terms with different n) . Each contribution is 
the product of two multipliers describing different optical 
processes. The second multiplier describes Rabi oscilla- 
tions with the frequencies v n = 2gy/n + 1, occurring in 
the QD-chain, while the first one shows that the region of 
Rabi oscillations moves therewith in space with velocity 
v = ^2a 2 k. One should note, that in the case of the exact 
synchronism v does not depend on n. That is why in this 
partial case the quantum statistics of light is not distorted 
with propagation of the wavepacket. In particular, the 
initial coherency of quantum light persists in time. It is 
immediately follows from Eq. (|34[) and normalization con- 
dition that the integral inversion w(t) corresponding to 
the inversion density (|45p oscillates in time similarly to 
the Rabi-oscillations in the single two-level system (see 
Eq.(6.2.21) in 1 at A = 0). 

Generally, each pair of the subpackets propagates due 
to the tunneling through its own potential barrier. The 
velocities of the subpackets motion are proportional to 
the barrier transparency (vf 2 ~ £1,2). If one of the bar- 
riers becomes absolutely opaque, the corresponding pair 
of subpackets does not move: vf 2 — ► at £i i2 — ► 0. 
The light propagation along the QD-chain is another nec- 
essary condition of Rabi subpackets motion: vf 2 — > 
at k — ► 0. 

For more detail analysis of Rabi wavepackets dynam- 
ics (in particular for taking into account the diffraction 
spreading) the integrals ([3"T)) - ([3"5)) need to be calculated 
numerically. The results of calculations is analysed in the 
following section. 



C. Coherent state: collapses and revivals picture 

Consider the case of an excited initial state with space 
distribution in form of single Gaussian beam: A n (x, 0) = 

c(n) exp(— x 2 /2cr 2 )/v / 7rtr 2 , B n+ x(x,0) = 0. Assume that 
initially light is prepared in coherent state, so pho- 
ton distribution is given by the Poisson law: c(n) = 

(n) n ^ 2 e~^ I 2 1 \[n\, where (n) defines the average pho- 
ton number. The spatial-temporal dynamics of inversion 
density for this case is depicted on Fig[4ji. As is seen from 
it, original Gaussian packet oscillates in time and moves 
along the chain. Oscillations collapse to zero quickly, 
but after a while they revive in another area of space. 
The phenomenon of Rabi oscillations collapses and re- 
vivals, caused by discreteness of photon distribution, is 
well studied^, but the spacing the collapse and subse- 
quent revival is qualitatively new effect. 

Although a variation of the inversion density, depicted 
in Figd^,, in arbitrary point of the space occupied by the 
Rabi wavepacket is not too large, the integral inversion, 
presented in FigHp oscillates between -1 and 1, thus in- 
dicating presence of strong light-QD coupling. 

Let us consider now the case of a mixed initial state 
(both A n (x,0) and B n+ x(x, 0) are nonzero). As it has 
been noted above, the synchronism condition can not 
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FIG. 4: Space-time distribution of the inversion density (a) 
and temporal dependence of the integral inversion (b) in the 
QD chain for a coherent initial state of light ((n) = 4). The 
initial state of QD-chain is a single Gaussian wavepacket. 
A n (x,0) = c(n)exp(-x 2 /2a 2 )/tfu^, B n+1 (x,0) = 0, £i = 
Wg, £ 2 = 7g, A = 2(£i - £ 2 ) + £ 2 a 2 fc 2 , fca = 0.5, a = 20a, 
A = 0.055. 



be fulfilled for both pair of subpackets (except the case 
when one of the pairs is motionless, see below). As a 
consequence, we have four subpackets (or three ones if 
A e // = for one of the pairs). If £1 = £2 = £ (poten- 
tial barriers with equal transparency) and A = 0, then 
both packet pairs move with pairwise equal velocities in 
opposite directions (vf = — V2 ). If the initial spatial 
distributions of probability amplitudes for ground and 
excited states are spaced, the subpackets will collide (see 
Figj5]). Initial Gaussian profiles of packets deform with 
time due to the difference in subpacket velocities. 

One should note that for the case of quantum exter- 
nal field in coherent state QD-chain cannot be saturated, 
i.e. the integral inversion could not be a constant (see 
Fig. [6|. This is due to the fact that A n (x,t) and 
B n+ i(x,t) are asymmetrical with respect to n and as 
a result E„ X!L I Ax (a, t)\ 2 dx - £ n \B n+1 (x,t)\ 2 dx 
is not a constant for any form of space dependence of 
A n (x,t), B n+ i{x,t). This asymmetry exists due to the 
presence of vacuum term in the sum and increases with 
decreasing of the average photon number (n) (compare 
FigGDand FigED. 

In other limiting case, for £2 = (fully opaque ground- 
state barrier), the velocities = and one pair of 
subpackets do not moves along the chain. If, moreover, 
A e ff(k/2) — 0, then = u£~ , and the synchronism con- 
dition is fulfilled for both packet pairs (see Fig[5]). 




20 150 



FIG. 5: Space-time distribution of the inversion in the QD 
chain for the case of two counterpropa gatin g wavepackets. 
A n {x,0) = c(n) exp[-(z - 3o-) 2 /2a 2 ]/ y^r^, B n+1 (x,0) = 
c(n + I) cxp[-(x + 3ct) 2 /2<t 2 ]/ v^Trcr 2 . The light is initially in 
coherent state ((n) = 5), A = , £1 = £2 = Wg, ka = 0.33. 
a) a = 10a, A = 0; b) a = 20a, A = O.O55. 



D. Fock qubit state and vacuum Rabi waves 

Let us now consider the interaction of the QD-chain 
with another important state of initial electromagnetic 
field, namely, the Fock qubit state, i.e. a superposition 
of two Fock states: |V/(0)) = Ci \N) +C 2 \N + 1), where 
N is an arbitrary fixed number, C\.i satisfy the normal- 
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FIG. 6: Temporal dependence of the integral inversion for 
a coherent initial state of light at the input parameters as 
follows: (n) = 5, a = 20a, A = 0.05g (solid line); (n) = 5, 
a = 10a, A = 0(dotted line) (n) = 0.5, a = 20a, A = 0(dashed 
line). In all cases, £1 = £2 = 10g, A = 0, ka — 0.33. 
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FIG. 7: The asymmetry of the inversion space-time distribu- 
tion for two counterpropagating wavepackets in case of the 
coherent state of light with small average number of the pho- 
tons (n) =0.5. A n (x,0) = c(n) exp[-(a;-3(T) 2 /2cr 2 ]/v / 47ra 2 , 
B n+1 {x,0) = c(n + l)exp[-(x + 3a)72a 2 ]/^I^, ft = ft = 
Wg, A = 0, ka = 0.33, a = 20a, A = 0. 




FIG. 8: Space-time distribution of the inversion in the QD 
chain for the case of the fully opaque ground-state barrier. 
A n (x,0) = c(n) exp[-(x - 3a) 2 /2a 2 ] / , B n+1 (a5,0) = 
c(n + 1) exp[-(x + 3a) 2 /2a 2 ]/^4^ 2 , ft = Wg, ft = 0, A = 
2(6 - 60 - £ia 2 fc 2 , = 0.5, a = 20a, A = 0.05g. The light 
is initially in the coherent state with (n) = 5 
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FIG. 9: Space-time distribution of the inversion (a) and tem- 
poral dependence of the integral inversion (b) in the QD 
chain for the field in the Fock qubit initial state. \ipf(0)) = 
1/V2|0> + 1/V2|1>, A (x,0) = Ai(as,0) = 
exp(-i 2 /2cr 2 ), A n (x,0) = for n > 2, B n (x,0) = for 
all n, Ci = 10<?, 6 = 7g, A = 2(^ - £ 2 ) + &a 2 k 2 , a = 20a, 
ka = 0.33, A = 0.05p. 



cited QD-chain is illustrated at Fig[9j This excitation 
couples single mode electromagnetic vacuum state with 
one-photon state only and can be imagined as a wave 
beam which is described by the monochromatic Rabi- 
frequency spectrum and continuous spatial spectrum at 
the same time. However, it should be noted that one 
need to consider an interaction of the QD-chain with a 
multi-mode field for full treatment of such effects^ . This 
problem is a subject for future considerations. 



ization conditions. As in the previous case, for initial 
spatial distributions the Gaussian beams will be used. 
The spatial-temporal dynamics of the inversion density 
for the initial state of field |V>/(0)) = 1/V2 |0) + 1/V2 |1) 
is shown in Fig. [9^. It can be seen from the figure that 
Rabi wavepacket does not collapse and revive, but oscil- 
lates in a complicated manner (see also integral inversion 
dynamics, FigJ5]). 

In the end of this section we consider the vacuum initial 
state of external electromagnetic field. It is well known 
that for this case the vacuum harmonic Rabi oscillations 
occuring due to spontaneous emission can take place 1 . 
The corresponding Rabi-wave packet in the initially ex- 



V. CLASSICAL LIGHT LIMIT 

The case of a coherent state of electromagnetic field 
with the large photonic number n is especially interest- 
ing for us. In this case one can neglect the quantum 
nature of the light and consider electromagnetic field as 
a classical one replacing the electric field operator with 
its expectation value E(x,t) — Re{<£ exp[i(kx — uit)]}. 
In doing this we neglect any changes of the field state. 
Omitting the local-field effects as before, we can write the 
Hamiltonian of the system " QD-chain - electromagnetic 
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FIG. 10: Space-time distribution of the inversion in the QD 
chain for the vacuum initial state of the field. Ao(x,0) = 
(1/v 7 ^) exp(-x- 2 /2cr 2 ) , A n (x,0) = for n > 1, B n {x,0) = 
for all n, £i = 6g, & = 5g, A = 2(£i + &a 2 k 2 , a = 10a, 
fca = 0.5, A = 0.05p. 



field" in the form H = Hq + Ht, where the term 



H f) = 



o \ - 



9 



+ ^i(nka—ut) _|_ jj 



(46) 

describes Rabi oscillations in non-interacting QDs and 
= fj,£/h is the Rabi frequency^, The interdot inter- 
action mechanism is independent on light properties, so 
the term Ht is defined as before by equation (J3J). 

The state vector of the system has the form: 



\*(t)) = J2(Mt)\ 



B P (t) 



(47) 



where A p (t), B p (t) are unknown functions. Now the 
equations of motion for them can now be obtained us- 
ing approach of Sec lII Bl Introducing the decay factor A 
into the nonstationary Schrodinger equation and making 
continuous limit transition in the QD-chain similar to the 
Sec lII CI we obtain the system of equations as follows: 



ft* = 



O f iUJ ° - X f , Ml *\ 

2i£ —a z - -I + —— n[kx - ut) 



(48) 



oo 

x J p(h,tMh,0)e^ hx+ ^^ + ^ h ^/ 2 Uh. 

— OO 

(49) 

where the values p(h,t), </5 ± (/i,t), iph{t), S±(h), Sl(h), 
A e ff(h) are determined, respectively by equalities (|26[) . 
p3]) . pOj) . (fT6j) . (fT7|) if one makes the substitution 



2g\/n + 1 — > r2^j (in doing that the photonic number de- 
pendence in p(h,t), ip , tph, ft disappears). The vector 
function 3?(/i, 0) is determined by the initial conditions 
similar to Eq. (f5Tj) . 
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FIG. 11: Space-time distribution of the inversion in the QD 
in the classical light chain for a single Gaussian wavepacket 
A(x,0) = (l/v^S 2 ) exp(-x 2 /2a 2 ), B(x,0) = 0, fi = 3fi fl , 
£2 = 0.9&, A = 2((, 1 -&)+&a 2 k 2 , ka = 0.33, 2n\/Q. R = 0.1. 
a) a = 5a; b) cr = 15a 



where *&(x,t) = 



Note, that spontaneous 



A(x,t) 
B(x,t) 

emission can be taken into account in the same way as 
in III Bt i.e. by describing it as interaction of QDs with 
photon bath. In this case the damping constant A in 
Eq. (|47|) will be the sum of two partial constants. The 
system (|47p can be considered as a particular case of the 
general system (jlip and can be solved in the same way. 
The solution reads 



A typical space-time distribution of the inversion den- 
sity w(x,t) = a[\A(x,t)\ 2 - \B(x,t)\ 2 } is shown in fig. [Til 
As is seen from the figure, the space-time dynamics of 
Rabi wavepacket is similar to the one for the case of vac- 
uum oscillations (but it is of substantially different physi- 
cal nature: electromagnetic field is in coherent state with 
photon number n — > oo and this state does not change 
when the system oscillates between ground and excited 
states in contradiction to the case of vacuum Rabi oscilla- 



tion). Since the photon distribution does not change, col- 
lapses and revivals are absent. However there is another 
phenomenon of similar nature. One can see in fig. I12[ 
where the plots of the integral inversion are presented, 
that the oscillations of this quantity at k ^ 0, £1 7^ £2 
decrease with time even at A = 0, whereas such a damp- 
ing is absent at k = 0, £1 = £2 and integral inversion 
oscillates harmonically in the range from -1 to 1 (dotted 
curve in fig. 1 1 2() . Such a behaviour is the consequence of 
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FIG. 13: Space-time distribution of the inversion density in 
the QD chain in the classical light for two counterpropagat- 
ing identical Gaussian wavepackets: A(x, 0) = exp[— (x — 
3a) 2 /2a 2 ]/^H^, B(x,0) = exp[-(x + 3a) 2 /2a 2 }/ Vi^ 2 , 
A = 0, ka = 0.33, a = 5a, £1 = £2 = 3fij{, ^ttX/Qr = 0.1. 



FIG. 12: Temporal dependence of the integral inversion in 
the classical light for the input parameters as follows: A = 0, 
k = 0, fi = £2 = 3Qr (dotted line); A = &a 2 k 2 , ka = 0.33, 
£1 = 6 = 3fi fl (solid line); A = 2(£i-,£2)-K2a 2 fc 2 , fca = 0.33, 
£1 = 3£Ir, £2 = 1-5£Ir (dashed line). In all cases a = 5a, 
A = 0. 

dependence of the effective detuning A e ff(h) (and there- 
fore the Rabi frequency f2(/i)) on h. But if k = and 
£1 = £2, then A e // = A and f2 does not depend on h. 
The damping rate is controlled by values of wavenumber 
k and coupling constants £1,2- If the initial state of the 
system is the excited state (B(x,0) = 0), the damping 
rate is predominantly determined by product £2^ (and 
damping is practically absent if £2^ = 0), for the ground 
initial state (A(x, 0) = 0) it is defined by £i/c. In the case 
of a mixed initial state all three quantity are important. 

Note, that more narrow wavepacket spreads faster, 
than the wider one (compare fig. [TTa and fig. [TTb). 

In the weak coupling limit the indicated dcphasing 
mechanism is analogous to the Landau damping in 
plasma. 

Interaction of two counterpropagating identical Gaus- 
sian Rabi wavepackets colliding at x — is shown in 
fig. [T31 Unlike the case of quantum field, the space-time 
distribution of the inversion density is fully symmetrical 
with respect to the line x — and integral inversion does 
not oscillate: this quantity equals zero for all t > and 
arbitrary values of Qr. 

VI. INFLUENCE OF THE LOCAL-FIELD 

EFFECTS 

Local fields arc depolarization fields arising inside QDs 
due to the dipole-dipole electron-hole interaction. As a 
result, the field acting on the e-h pair differs from mean 
field^. The local-field action is described by Eq.(0 and 



leads to nonlinear terms in the equations of motion (last 
terms in Eqs. ([6]) -10). The presence of nonlinearity re- 
sults in the qualitatively new features of Rabi oscillations 
in single QD o 17 i 18 ' 19 . The influence of the local-fields on 
Rabi waves is the subject for future investigations. In 
this section we consider some aspects of this influence 
only and discuss some sufficient conditions of the local- 
field negligibility. 

Let us consider finite (of the length L) QD-chain inter- 
acting with classical light. The relaxation processes are 
neglected (i.e. A = 0). The equations of motion follow 
from Eqs.([6])-(0 in the same way as Eq. fTTj) reduce to 
Eq. pE|) and have the form 

d t A = - 40 A + l ^Be^-^ + 

i£,a 2 d 2 x A-iAoj\B\ 2 A, (50) 

d t B = -(u> + 4i)B + t lh Ae -i(kx-wt) + 

i£,a 2 d 2 x B -iAu\A\ 2 B . (51) 

We use the periodic boundary conditions (Born-von 
Karman conditions): 

A(L/2,t) = A(-L/2,t), B(L/2,t) = B(-L/2,t). (52) 

Suppose that electromagnetic field are also satisfy peri- 
odic conditions, i.e. kL = 2 / nm 1 where m is integer num- 
ber. Then the problem under consideration has a clear 
physical interpretation: we are dealing with ringed QD- 
chain of the radius R = L/2n. QD-chain interacts with 
electromagnetic mode of whispering gallery type, which 
is excited in microcavity with rotation axis ( e lkx — > e LmLp , 
x — > Rip, ip is the azimuthal coordinate). 

Let us seek partial solution of the sys- 
tem ([50)1 - (f5T|) in the form of traveling wave: 
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A(x,t) — uoe 

i(h-k/2)x p -i(v-u)/2)t 



i(h+k/2)x -i(v+u>/2)t 



B(x,t) = 

Wo g-v» "\- "i-i-^ w here uo, Vq are constant 

amplitudes satisfying the normalization condition for 
wavefunction, h and v are real values. Ssubstituting A 
and B into Eqs. (f50]) - (f51~j) we obtain 



21 

Aw|i> | J wo + — uq = 0, 



9. 



\v + 02 - Awltiol^ I "0 + ^Tr^O = 



(53) 
(54) 



where </>ij(h) = 2£ 1>2 T A/2 - £i, 2 a 2 (/i ± fc/2) 2 . From 
boundary conditions (|52[) one can obtain the quantization 
condition for h: h — im/L, n is integer. Expressing vq 
from Eq. ([54| and substituting into Eq. (f53|) . we obtain 



i(h) 



Sl 2 a (v + fa(h)) 



A{v + fo{h) - Auj\utff 
The wavefunction normalization condition results in 



0. 



(55) 



1 



n 2 



4\u + fa(h) - Aui\u G \ 



212 



"0 = T 



(56) 



Eqs. ([55l) - ([56|) form the closed system of equations with 
respect to v and |ito| 2 , which defines the dispersion law of 
Rabi waves v — v(h) — v(n) with regard to the local-field 
effects. It should be noted that only solutions with real- 
valued v have physical meaning (otherwise the system 
([53")) - ([54")) contradicts to the initial equations ([5D)) - (|5T|) V 
For L — ► oo, Au — * the system (l5"5 j) - (|5o]) reduces to 
the dispersion equation^ for Rabi waves in the infinite 
QD-chain interacting with classical light. 

The dispersion characteristics of Rabi waves for two es- 
sentially different values of L but identical Aui and other 
parameters are shown in Fig ll4l The curves depicted 
in FigfTib corresponds to the QD-chain of large radius. 
The nonlinear effects in this case are negligibly small. 
Dispersion curves do not feel nonlinear effects and are 
practically identical to obtained in 45 . For QD-chain of 
small radius the situation becomes principally different 
due to nonlinear effects. The closed loop appears in the 
vicinity of any cross point and the number of the modes 
increases to four (FigfMbV The size of this loop decreases 
with decreasing the ratio Au>/£Ir. 

The detailed structure of closed loop is shown in 
FigfHb. When h (and, correspondingly, n) changes adi- 
abatically slow the motion along the loop takes place (it 
is shown by arrows in Fig[T4"b'). The transitions from one 
branch to another are followed by simultaneous changes 
of Rabi frequency and Rabi wave group velocity. The 
trajectories of motion in different directions proves to be 
nonidentical. Thus, the loop has an hysteresis character. 
Interaction of the local-fields and Rabi waves leads to ef- 
fect of bistability of quantum states in the QD-chain. The 
detailed analysis of different dispersion branches stability 
is the subject for future investigations. 

Let us now discuss the conditions of the local field 
effects neglecting in equations Note, that for 




FIG. 14: Dispersion curves for Rabi waves subject to local- 
field effects, a) L = 200a, b), c) L = 20a. In all cases 
£l = 6 = 3fi fl , ka = 0.25tt, A = 0, Auj = 100fi H . 



Ci.2 — > the system (HJ)-© transforms to equations 
(24) ir>i2, where Rabi oscillations in single QD are de- 
scribed with account of local field effects. An analy- 
sis of equations (24) ir>i^ has shown that local field ef- 
fects are negligibly small for the case of sufficiently large 
Rabi frequencies, namely, f2/ n \ > Auj, where fi/ n \ = 
-y/A 2 + Ag 2 ((n) + 1), (n) is average photon number. If 
this condition is fulfilled for coherent state, the con- 
ventional picture of Rabi oscillations (the alternation 
of collapses and revivals) takes place. In the opposite 
case qualitatively new regime of oscillations arises. This 
regime is due to the local field effects and is described 
in detail in 1 ^. In this work we are interested only in 



first regime. In our case the role of 0/ n \ is played by 
frequencies vi^in^h), which are determined by Eq. (|15|) . 
Therefore it is reasonable to suggest that for traveling 
Rabi waves the local fields effects are negligible under 
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the condition 



^1,2 (n, h) > Aiv, 



(57) 



where Aoj is determined by Eq.®. 

In the case of Rabi-wavepackets general line of reason- 
ing is the same as for Eq. (|57f and leads to the following 

condition of the local-field negligibility: vi t 2 ((n) , h^l) > 
Auj. 



In this section we will calculate exciton-exciton and 
exciton-photon multi-time correlation functions. For this 
purpose it is convenient to pass to the Heisenberg picture 
and to use the discrete model of QD-chain (transition to 
the continuous limit will be made on the final stage). 
Thus, we should look for a time dependence of operators 
(7^, i.e. we shall seen to find the evolution operator. It 



can be obtained from 
and has the form 



(for details see Appendix [2J 



VII. CORRELATION FUNCTIONS 




■ 

I 



FIG. 15: The space-time dependence of normalized corre- 
lation function (x, 0, t, 0) for light in Fock state with 
n — 5 and the input parameters as follows: £i = £2 = Wg, 



A 



-i(x,0) = ( 



5a, A„(a;,0) = exp(-x A /2a 2 )/^ 



7T / a 

^ ) = |^EE / {<Pn(h,t)e< t \a p ,n)(a q ,n\ + 



— re/a 

ip n (h,t) e 4 ^* \ap,n) {b q ,n + l\ + e iS ^ t |6 p ,n+l) (a q ,n\ 



+<p+(h, t)e iS ^ \bp, n + 1) (b q , n+l\\ e th ^ a dh. 



(58) 



Here ip^(h,t), ip n (h,t), £l n (h) are given by formulas 
([2"5]). pbj) . (fTf]|) . Quantities Sf, A eff {h) are defined 
by expressions ([30]) . (fT7|) if one substitutes #1,2 — > 
2£i,2 cos[(ft ± /c/2)a]. Note, that if one makes approxi- 
mation cos[(/i ± fc/2)a] = 1 - (h ± k/2) 2 a 2 /2 (this is the 
usability condition for continuous limit), the quantities 
5^, A e ff(h) coincide with expressions (f3TJ)) . (jTTj) . 



From expression (|58p we obtain by using the formula 
&p(t) = ^ + (^ 0)(7^ : L7(t, 0) the exciton-exciton correla- 
tion functions G p ]l(t, t') = (<7+ (t)a q (t')) as follows: 



G M^) = ^iEE Iff ^'^V^-^e^-^e-^*-^ [ VJ +(ft / ,t)^-(A',0-^(^,*)V' n (^,0] 

n Z,r -7r/a 

1 ( /l > *)Pn+l (f> *'K,n+l (°H,n+l (0) + ^n+l(h, t)tp n+X (g, t')«*,„+2 (0H,n+2 (0) (59) 

+ ^n+iC 1 ' *)^n+i (5, *') u f, n+i (0H,n+2 (0) - V'n+i (h, t)<p~ +1 (g, t')v* ln+2 (0)u r , n+1 (0)]dti dhdg , 

I 



where u p , n (t) = A p , n {t)e-^ a -^' 2 , u P ,„+i(t) = 
Sp,„ + i(i)e i (' c?>a - Wo *)/ 2 . The continuous limit transition 

can be made by means of the following substitutions: 

00 00 

pa .x, <?a x\ ~* S / •■■ rf2 -E~ > i / •■■ cfe/ - 

/ — 00 r —00 

One can see, that in general case G p ]l(t, t') can not be 

represented as G p ^ 9 (i, t') or Gp) q (t — t'). It indicates the 
absence both spatial and time homogeneity. 

Let us consider case of -B P) „ + i(0) = and t' — 0. Then 
ip n (h,t') = 0, <Pn(h,t') — 1 and taking into account that 

J-tt/o el9 ^ r ^ a dg — 2Tr5 qr /a, we obtain rather simple ex- 



pression for correlation function: 



°) = E E JJ V^'e** 1 -^ 

" I -iv/a 

ip+(h',t)cp+ +1 (h,t)ul n+1 (0)u q , n+1 (0)dh'dh, (60) 
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or, in continuous limit, 

G (1 \x,x' ,i,0) = 



a 2 e~ lAt 



ir/a 



/ JJ" e ih'(x-x') e ih(z-x) 
™ -oo -t/o 



8tt 3 

t)v>+ +1 (fc, tK +1 (*, 0)u n+ i(a:', 0)dh'dhdz. (61) 



Plot of normalized correlation function 

(!)/■ / + _ G^ (x,x' ,t,t') 



gW(x > x',t > t') = 



v /G( 1 )(x,a;,0,0)G( 1 )(x',x',0,0) 



for x' = 0, f' = is presented on Fig|15] The external 
field is supposed to be in a single-photon initial state and 
the electron in QD-chain is in the excited initial state 
(B n+1 (x,0) =0). 



The correlation function of the polarization operator 
a~p(i) with the operator of the external field E~(t) — 



(62) £de^ kpa (exciton-photon correlators) appears as 



£ 

^(^(^i^EE If dhdh'{^ +1 {h\t)V^[^ n {h,t)ul n ^ (63) 

n (,r -7r/a 

-^„ +1 (/^)V^T2[^(MK„+i(0H„(0) + ^+(/ l ,t)< n+1 (Q) Ui , n+1 (0)]}e i ^(' l - ft ')e i ( /tV - w ) a e i ( fc ^- Wot )e^(*- t '). 

I 



7T / a 



VIII. CONCLUSION 

The main results, obtained in this work, can be sum- 
marized as follows: 

i) In QD-chain with tunneling coupling in strong 
electron-photon coupling regime the space propagation 
of Rabi oscillations (Rabi waves) takes place. For prop- 
agation of Rabi waves the wave vector of the photon 
mode must have a nonzero component along the chain. 
Characteristics of the Rabi waves depend strongly on re- 
lations between parameter of electron-photon coupling, 
frequency deviation and transparency factors of tunnel- 
ing barriers for both of levels. 

ii) Traveling Rabi waves are quantum states of the QD- 
chain dressed by radiation. The qualitative distinction 
of this states from states of single dressed atorn^ is the 
space-time modulation of dressing parameter by travel- 
ing wave law. Traveling Rabi waves can be interpreted 
as entangled states of e-h pair and photons. These states 
are characterized by the dependence of energy on quasi- 
momentum and can be treated as new type of quasipar- 
ticles (in this paper they are called rabitons). 

iii) There are two traveling Rabi modes with differ- 
ent frequencies of Rabi oscillations at a given value of 
wavenumber. The range of Rabi oscillations frequencies 
is limited by the critical value, which is different for both 
modes. The QD-chain is opaque for Rabi modes with fre- 
quencies lesser than critical one. The critical frequencies 
as well as dispersion characteristics depend on number of 
photons in the mode. 

iv) In general case the propagation of Rabi wavepack- 
ets occurs in the form of four partial packets, which are 
characterized by different amplitudes and velocities of 
motion. Rabi wavepackets transfer energy, inversion, 
quasimomentum, electron-electron and electron-photon 
quantum correlations along the chain. Particularly, in 



the case of QD-chain interacting with quantum light in 
the coherent state the known collapses-revivals picture^ is 
drastically modified due to propagation effect: collapses 
and revivals take place in different points of space. 

v) One pair of Rabi wavepackets exists due to the 
ground state contribution in initial state of the system, 
another one appears due to the contribution of excited 
state (generally, initial state is arbitrary superposition of 
ground and excited states). Each pair of packets moves 
through its own potential barrier. If one of the barriers is 
completely opaque, then corresponding pair of wavepack- 
ets does not move and does not spread. 

vi) Rabi wavepackets propagation along the QD-chain 
is followed by the transformation of quantum light statis- 
tics (for example, in initially coherent light incoherent 
component appears). If specific conditions are fulfilled 
then incoherent part can be decreased to zero. 

vii) In particular cases a number of Rabi wavepackets 
can be diminished. There are two mechanisms of such 
diminution: tending to zero of partial wavepacket ampli- 
tude and merging of the packets due to their velocities 
equality. 

Rabi waves can take place in a number of other dis- 
tributed systems strongly coupled with electromagnetic 
field. 

The example are superconducting circuits based on 
Josephson junctions, which are currently the most exper- 
imentally advanced solid-state quintal. It is evident, for 
example, that the chain of qubits placed inside a high-Q 
transmission-line resonator due to the qubit-qubit capac- 
itance coupling will support the Rabi waves propagation 
similar to described in this article. 

Rabi waves effect can be practically used in nanoelec- 
tronics, quantum computing, quantum informatics. 
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APPENDIX A: EVOLUTION OPERATOR 

In this appendix, the evolution operator for a dis- 
crete QD-chain fEq. (j5"8|) ) is obtained from Eq.© Let 
us introduce new variables u p .„ = A p , n e~ l ^ kpa ~^ ot ^ 2 , 
v p ^ n+ i = B p , n+ ie i{ - kpa -^° t)/2 . It leads to the system 



a ■ t ( —ika/2 , ika/2 



- ig\/n + lv p .n+ie iAt , (Al) 
d t v p , n+1 = i& (v p - 1 , n+1 e lka/2 + v p+hn+1 e- lka / 2 ^ 

- igVnTTu p . n e~ iAt . (A2) 



Coefficients in this system do not depend from QD- 
numbcr p. It allows us to seek the solution in the 
form Up±i,„ = u p , n e ±lha , v p ±i,n+i = w p ,„ + ie ±l ' ia , ha € 
[— 7r, 7r]. Then the equations for different QDs become 
independent and may be written as 

dtUp n = 2t£i cos[(/i + k/2)a]u p n — ig\J n + lv p „+ie lAt , 

(A3) 



d t Vp.„ + i = 2i£ 2 cos[(/i - k/2)a]v p , n+1 - ig^n + lu p . n e 

(A4) 



-iAt 



Expressing v p . n+ i from (|A3|) and substituting into (IA4|1 . 
one obtains the set of independent second-order ordinary 
differential equations for the variables u Pi n 



d V'p.n . c dllp n 
lOh- ' 



dt 2 



d.t 



(A5) 



where S h = A + 2£i cos[(/i + k/2)a] + 2£ 2 cos[(/i - k/2)a]. 
The solution of (IA5D reads 



up n (t) = c™(o)e j[ ^+°" ( ' i)It/2 + c^oyt^-^* 1 ^/ 2 

(A6) 

with the constants C™ 2 (0) to be determined from ini- 
tial conditions. Writing the latter in the form Up,„(0) = 
a n e lhpa : v p , n+ i(0) = b n+1 e lhpa , we rewrite (JXg) as 



2 1 



and 



u ihpa 

(A7) 



^ihpa 



(A8) 



A, 



gv^ + 1 , 

fi»(fc) 



n n (h) + A h 

2Q n (h) 
Cl n (h) - A h _ gVn + 1 

2Q n (/*) °" n n (h) 0,1+1 ' 



n n (h) 



2n n (h) 



b n -\ 



gy/n + l n n (h) + A h 

-0,n H — — — On+1- 



n n (h) 



2fl n (h) 



(A9) 
(A10) 
(All) 
(A12) 



Basis wavefunction has the form 



a p , n) + t>p, n +i(X) IV + 1)) ■ 



n p 

(A13) 

We now can represent the required wavefunction as 
Fourier integral with basis functions ^ft(i): 

7r / a 

!*(*))= / M(fe)|* h (t))dfe, (A14) 

— 7r/a 

where M(h) is unknown weighting function. To deter- 
mine it let us take into account, that 



7r /a 

l*(°))=E / ^)K|« P ,™)+&n + i|6 P ,«+l)]e ttp "rfft, 

(A15) 



— 7T /a 



where 



7T / a 

(a p ,n |¥(0)) = a„ / M(h)e lhpa dh (A16) 

— 7r/a 

and the same for (6 p ,n + 1 1^(0)). Then employing in- 
verse Fourier transform we have 



„M(h) = ^ J2 (« 9 ^l*(0))e-^ a (A17) 



q— — oo 



and analogously for b n+ \. 

Substituting (|aT7| into ([AT4]) . we obtain |*(t)) = 
U(t,0) |*(0)), with L> given by the Eq.(|5SJ). These cal- 
culations allow one to pass from Schrodinger picture to 
the Heisenberg one. They are used in the calculations of 
different types of correlators. 
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